How to Make the Quantum Adiabatic Algorithm Fail 



Edward Farhi/'Q Jeffrey Goldstone, 1 Sam Gutmann, 2 and Daniel Nagaj 1 

1 Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts, 02139 
2 Department of Mathematics, Northeastern University, Boston, Massachusetts, 02115 

(Dated: Dec 19, 2005) 

The quantum adiabatic algorithm is a Hamiltonian based quantum algorithm designed to 
hnd the minimum of a classical cost function whose domain has size N . We show that poor 
choices for the Hamiltonian can guarantee that the algorithm will not hnd the minimum if 
the run time grows more slowly than >/~N. These poor choices are nonlocal and wash out 
any structure in the cost function to be minimized and the best that can be hoped for is 
Grover speedup. These failures teli us what not to do when designing quantum adiabatic 
algorithms. 



I. INTRODUCTION 

The quantum adiabatic algorithm was introduced 0] as a quantum algorithm for finding the 
minimum of a classical cost function h(z), where z = 0, . . . , N — 1. This cost function is used to 
define a quantum Hamiltonian diagonal in the z basis: 

AT-l 

H P = £ h(z)\z). (1) 

2=0 

The goal is now to find the ground state of Hp. To this end a "beginning" Hamiltonian Hp is 
introduced with a known and easy to construct ground state \gp). The quantum computer is a 
system governed by the time dependent Hamiltonian 

H(t) = (1 - t/T)H B + (t/T)H P , (2) 

where T controls the rate of change of H(t). Note that H(0) = Hp and H(T) = Hp. The state of 
the system obeys the Schròdinger equation, 

ij t m)) = H(t)\mh (3) 

where we choose 



1^(0)) = \9b) 

and run the algorithm for time T. By the adiabatic theorem, if T is large enough then \tp(T)) will 
have a large component in the ground state subspace of Hp. (Note we are not bothering to state 
the necessary condition on the lack of degeneracy of the spectrum of H{t) for < t < T, since it 
will not play a role in the results we establish in this paper.) A measurement of z can then be used 
to find the minimum of h(z). The algorithm is useful if the required run time T is not too large as 
a function of N. 

There is hope that there may be combinatorial search problems, defined on n bits so that N = 
2 n , where for certain "interesting" subsets of the instances the run time T grows subexponentially 
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in n. A positive result of this kind would greatly expand the known power of quantum computers. 
At the same time it is worthwhile to understand the circumstances under which the algorithm is 
doomed to fail. 

In this paper we prové some general results which show that with certain choices of Hp or Hp 
the algorithm will not succeed if T is o(y/N), that is T/-/N — > as N — > oo, so that improvement 
beyond Grover speedup is impossible. We view these failures as due to poor choices for Hp and 
Hp, which teach us what not to do when looking for good algorithms. We guarantee failure by 
removing any structure which might exist in h(z) from either Hg or Hp. By structure we mean 
that z is written as a bit string and both Hp and Hp are sums of terms involving only a few of 
the corresponding qubits. 

In Section II we show that regardless of the form of h(z) if Hp is a one dimensional projector 
onto the uniform superposition of all the basis states \z), then the quantum adiabatic algorithm 
fails. Here all the \z) states are treated identically by Hp so any structure contained in h(z) is 
lost in Hp. In Section III we consider a scrambled Hp that we get by replacing the cost function 
h(z) by h(jr(z)) where tt is a permutation of to N — 1. Here the vàlues of h{z) and h{-ïï{z)) are 
the same but the relationship between input and output is scrambled by the permutation. This 
effectively destroys any structure in h{z) and typically results in algorithmic failure. 

The quantum adiabatic algorithm is a special case of Hamiltonian based continuous time quan- 
tum algorithms, where the quantum state obeys © and the algorithm consists of specifying H(t), 
the initial state \ip(0)), a run time T and the operators to be measured at the end of the run. In 
the Hamiltonian language, the Grover problem can be recast as the problem of finding the ground 
state of 

H w = E(I-\w)(w\), (4) 

where w lies between and N — 1. The algorithm designer can apply H w , but in this oracular 
setting, w is not known. In reference [j| the following result was proved. Let 

H(t) = H D (t)+H w , (5) 

where Hp is any time dependent "driver" Hamiltonian independent of w. Assume also that the 
initial state \i/j(0)) is independent of w. For each w we want the algorithm to be successful, that is 
\ip(T)) = \w). It then follows that 




The proof of this result is a continuous-time version of the BBBV oracular proof [3]. Our proof 
techniques in this paper are similar to the methods used to prové the result just stated. 



II. GENERAL SEARCH STARTING WITH A ONE-DIMENSIONAL PROJECTOR 

In this section we consider a completely general cost function h(z) with z = 0, . . . , N — 1. The 
goal is to use the quantum adiabatic algorithm to find the ground state of Hp given by Q with 
H(t) given by ©. Let 

1 

\s) = -7= V \z) (7) 
y/N ^ 

v 2 = 

be the uniform superposition over all possible vàlues z. If we pick 

H B = E(I-\s)(s\) (8) 
and |^(0)) = \s), then the adiabatic algorithm fails in the following sense: 
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Theorem 1. Let Hp be diagonal in the z basis with a ground state subspace of dimension k. Let 

H(t) = (1 - t/T)E (I - \s)(s\) + (t/T)H P . 

Let P be the projector onto the ground state subspace of Hp and letb>0 be the success probabilíty, 
that is, b = (V>(T)|P|^(T)). Then 

r >i, ÍE ^ 



E\ k E 

Proof. Keeping Hp fixed, we introduce N — 1 additional beginning Hamiltonians as follows. For 
x = 0, . . . , N — 1 let V x be a unitary operator diagonal in the z basis with 

(z\V x \z) = e 2nizx / N 

and let 

V x \s) = ^=J2e^ N \z) 



N 

z=0 

so that the {\x)} form an orthonormal basis. Note also that 

\x = 0) = \s). 

We now define 

H x (t) = (1 - t/T)E(l - \x)(x\) + (t/T)H P , 

with corresponding evolution operator U x (t2,h). Note that H(t) above is Ho(t) with the corre- 
sponding evolution operator Uq. For each x we evolve with H x {t) from the ground state of H x (0), 
which is \x). Note that H x = V x H V x t and U x = V x U v£. Let \ f x ) = U x (T,0)\x). For each x the 
success probability is (f x \P\f x ), which is equal to b since P commutes with V x . The key point is 
that if we run the Hamiltonian evolution with H x backwards in time, we would then be finding x, 
that is, solving the Grover problem. However, this should not be possible unless the run time T is 
of order y/N. 

Let Ur be the evolution operator corresponding to an x-independent reference Hamiltonian 

H R (t) = (1 - t/T)E + (t/T)H P . 

Let \g x ) = ~^P\fx) be the normalized component of \f x ) in the ground state subspace of Hp. We 
consider the difference in backward evolution from \g x ) with Hamiltonians H x and Hp, and sum 
on x, 

S(t) = Y,\\u j x (T,t)\g x )-Uj i (T,t)\g x ) 

X 

Clearly S(T) = 0, and 

s (°) = E|^( T ' )i^)-^( T ' )i^)„ 

X 

Now \g x ) = Vb\f x ) + ^/T^b\f x l ) where |#) is orthogonal to \f x ). Since U x (T,0)\f x ) = \x) we have 

S(0) = \\y/b\x) + VT^blx^ - \i x ) 



where for each x, Ix -1- ) and \i x ) are normalized states with \x ± ) orthogonal to \x). Since Hp 
commutes with Hp, \i x ) = UÍ(T,0)\g x ) is an element of the fe-dimensional ground state subspace 
of Hp. We have 

5(0) = 2N - [Vb(x\i x ) + Vl - b{x^\i x ) + c.c. 

X 

> 2N - 2VbJ2\(x\i x ) - 2NV1 - b. 

X 

Choosing a basis for the k dimensional ground state subspace of Hp and writing \i x ) = 

a x i\Gi) H h a xk \G k ) gives 



y^ K^Ka;) < ^ \a x j\ ■ \(x\Gj)\ 



^ <u2\ a xj\ 2 Y.\( x '\ G i 



'Nk. 



Thus 



5(0) > 2N(1 - VT^b) - 2VbVM. 
We will use the Schròdinger equation to find the time derivative of S(t): 
d 



dí 



5(í) 



Y,^ [(9*\Ux(T,t)tf R (T,t)\g x ) + c.c. 

X 

iJ2(sMT,t)[H x (t) - H R (t)]U R (T,t)\g x ) + c.c. 

X 

2 Im (1 " t/T)E(g x \ U x (T, í) \x) (x\ U R (T, t) \g x ) . 



Now 



dí 



S(t) < 2E(l-t/T)Y\(9x\U x (T,t)\x)(x\U R (T,t)\g x ) 

X 

< 2E{l-t/T)Y J \{AU R (T,t)\g x 

X 

Using the same technique as in (jüj), we obtain 

< 2.67(1 - í/T)VM. 



Therefore 



Jo 



dí < ETVNk. 



Now 5(0) < 5(T) + J T \&S(t)\ dí and S(T) = so 



(9) 



(10) 



5(0) < ETVNk. 
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Combining this with (jlU|) gives 

ET\fWk > 2N(1 - y/l-b) - 2Vb\fWk, 
which implies what we wanted to prové: 




□ 

How do we intèrpret Theorem 1? The goal is to find the minimum of the cost function h(z) 
using the quantum adiabatic algorithm. It is natural to pick for Hb a Hamiltonian whose ground 
state is \s), the uniform superposition of all \z) states. However if we pick Hb to be the one 
dimensional projector E(I — \s)(s\) the algorithm will not find the ground state if T/yN goes to 
as N goes to infinity. The problem is that Hb has no structure and makes no reference to h(z). 
Our hope is that the algorithm might be useful for interesting computational problems if Hb has 
structure that reflects the form of h(z). 

Note that Theorem 1 explains the algorithmic failure discovered by Znidaric and Horvat 0] for 
a particular set of h(z). 

For a simple but convincing example of the importance of the choice of Hb, suppose we take a 
decoupled n bit problem which consists of n clauses each acting on one bit, say for each bit j 



hj(z) 



iízj = 0, 

1 ÏÍZj = 1, 



so 



h(z) = z\ + Z2 H h z n . (11) 

Let us pick a beginning Hamiltonian reflecting the bit structure of the problem, 

n 1 

H B = Y,2( l - a * ] )- (12) 

3=1 

The ground state of Hb is |s), The quantum adiabatic algorithm acts on each bit independently, 
producing a success probability of 

p=(l- q (T)) n , 

where q(T) — ► as T — > oo is the transition probability between the ground state and the excited 
state of a single qubit. As long as nq(T) — > const. we have a constant probability of success. 
This can be achieved for T of order ^/n, because for a two level system with a nonzero gap, 
the probability of a transition is q(T) = 0(T~ 2 ). (For details, see Appendix A.) However, from 
Theorem 1 we see that a poor choice of Hb would make the quantum adiabatic algorithm fail on 
this simple decoupled n bit problem by destroying the bit structure. 

Next, suppose the satisfiability problem we are trying to solve has clauses involving say 3 bits. 
If clause c involves bits i c , j c and k c we may define the clause cost function 

, . _ í if z ic , Zj c , z kc satisfy clause c, 
\ 1 otherwise. 
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□ projector 
o bit structure 



10 11 12 13 14 15 

n bits 



FIG. 1: Median required run time T versus bit number. At each bit number there are 50 random instances 
of Exact Cover with a single satisfying assignment. We choose the required run time to be the value of 
T for which quantum adiabatic algorithm has success probability between 0.2 and 0.21. For the projector 
beginning Hamiltonian we use (JSJ) with E = n/2. The plot is log-linear. The error bars show the 95% 
confidence interval for the true medians. 



The total cost function is then 



h{z)=Y,h c {z). 



To get Hb to reflect the bit and clause structure we may pick 



Hb,, 



( 1 _ ff W) + (l_ a ü.J) + (i_ .W) 



with 



Hb = Hb,, 



(13) 



In this case the ground state of Hg is again \s). With this setup, Theorem 1 does not apply. 

We did a numerical study of a particular satisfiability problem, Exact Cover. For this problem 
if clause c involves bits i c , j c and k c , the cost function is 



h ( z \ = j if z ic + % + z k c = 1, 
1 1 otherwise. 



Some data is presented in FIG. 1. Here we see that with a structured beginning Hamiltonian the 
required run times are substantially lower than with the projector Hb- 



7 



III. SEARCH WITH A SCRAMBLED PROBLEM HAMILTONIAN 



In the previous section we showed that removing all structure from Hp dooms the quantum 
adiabatic algorthm to failure. In this section we remove structure from the problem to be solved 
(Hp) and show that this leads to algorithmic failure. Let h{z) be a cost function whose minumum 
we seek. Let tt be a permutation of 0, 1, . . . , N — 1 and let 

hW(z) = Ii(tt~ 1 {z)) . 

We will show that no continuous time quantum algorithm (of a very general form) can find the 
mimmum of fcM for even a small fraction of all tt if T is o(y/N). Classically, this problem takes 
order N calls to an oracle. 

Without loss of generality let h(0) = 0, and h(l),h(2), . . . ,h(N — 1) all be positive. For any 
permutation tt of 0, 1, . . . , N — 1 we define a problem Hamiltonian Hp n , diagonal in the z basis, as 

N-l N-l 

Hp,, = h M (z)\z)(z\ = h(z)\n(z)){K(z)\. 

2 = 2=0 

Now consider the Hamiltonian 



H n (t) = H D (t) + c(t)H Pi7T (14) 

for an arbitrary 7r-independent driving Hamiltonian Hp>(t) with |c(í)| < 1 for all t. Using this 
composite Hamiltonian, we evolve the 7r-independent starting state \ip(0)) for time T, reaching the 
state \^ n (T)). This setup is more general than the quantum adiabatic algorithm since we do not 
require Hp,{t) or c(í) to be slowly varying. Success is achieved if the overlap of \tp w (T)) with 1 7r(0)) 
is large. 

We first show 



Lemma 1. 

\\\MT)) ~ WAT))f < 4h*TN\VW=ï, (15) 

7r,7r' 

where the sum is over all paírs of permutations tt, tt' that differ by a single transposítíon involving 
vr(0), and h* = ^J2Hz) 2 /{N-l). 

Proof. For two different permutations tt and tt' let |^/v(í)) be the state obtained by evolving from 
\ip(0)) with H n and let |?/v(í)) be the state obtained by evolving from IV'(O)) with H^/. 
Now 



d_ 

dí 



\Mt))-\Mt)) 



dt 



= i(ip n (t)\(H v {t) ~ Hv>(t))\^,(t)) + C.C. 



< 2 



<^ w (í)|(^(í)-^(í))|^(í)) 



We now consider the case when tt and tt' differ by a single transposition involving 7r(0). Specifically, 
tt' = tt o (a <-> 0) for some a. Now if tt(0) = i and tt(o) = j, we have tt'(0) = j and tt'{o) = i. 
Therefore, since /i(0) = 0, 



H P , n - Hpy = c(t)h(a) - = c(t)h(a) (|vr(a))(^(a)| - |7r'(a))(vr / (a)|) , 
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so we can write 
d 



dt 



E|llM*)>-hM*)> < 2\c(t)\J2Ka) <^(t)| (|7r(a)><7r(a)| - \Aa))(A")\) llM*)> 

tt,it' n,ir' 



This further simplifies to 



^£||hM*)> " hM*)>f < 2^>(a) (|(^(í)|7r(a))| + |<7r'(a)|^(í)>|) 

7r,7r' 7r,7r' 

= 2^^/ l (a)|(^(í)|7r(a))|+2^^%)|<7r»K,(í)> 

= 4^^%)|(^(í)|7r(a)) 

= 4^^%)|(^(í)|7r(a)) 



< 4^ /^/i(a) 2 = 4/i*iV!VÏV^Ï. 
7r y a 



where we used the Cauchy-Schwartz inequality to obtain the last line. Integrating this inequality 
for time T, we obtain the result we wanted to prové, 



]T \\\MT)) - MAT)) < 4h*TNlVN=ï, 



where the sum is over tt and w' differing by a single transposition involving 7r(0). 
Next we establish 



□ 



Lemma 2. Suppose |1), |2), \L) are orthonormal vectors and j(V'iK) 
where i = 1, . . . , L. Then for any normalízed \ip), 



^L·ii-lf) >bL -2VT. 



> 6 /or normalízed vectors 



(16) 



Proof. Write 



E||i^)-i^|| ^ £|<«-<#>| 

i i 



and use the Cauchy-Schwartz inequality to obtain 

2 



> bL-2\fl. 



□ 
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We are now ready to state the main result of this section. 

Theorem 2. Suppose that a continuous time algorithm of the form (|14JI succeeds with probabilíty 

2 



at least b, i.e. 



> b, for a set of eN\ permutations. Then 



e 2 b , 

T > ——y/N-1 - 
~ I6h* 4h* 



(17) 



Proof. For any permutation ir, there are N—l permutations 7r' a obtained from ir by first transposing 
and a. For each ir let be the subset of those N — l permutations on which the algorithm 
succeeds with probability at least b. Any such permutation appears in exactly N — 1 of the sets 
5^ so we have 

^2\S^\ = (N-l)eNl 
Let M be the number of sets S-^ with \S n \ > è(N — 1). Now 

^ ] | Sn | = ^ ] |«5tt| + ^ ] 

|5.[>f(AT-l) |5W|<5(JV-1) 

^1^1 < M(N-l) + (Nl-M)^{N-l), 



(iV-l)eiV! < M(iV- 1) + M-(JV- 1), 

so M > |iV!, i.e. at least |iV! of the sets «S^ must contain at least |(iV— 1) permutations on which 
the algorithm succeeds with probability at least b. For the corresponding 7r, we have 



£ IWCO) - |^(T)) > ò|(iV - 1) - 2 A /|(iV 



1). 



by Lemma 2. (Note that the algorithm is not assumed to succeed with probability b on ir.) Since 



there are at least |iV! such n, 



£ ||^(T)) - |^(r))f > \N\ [b e -(N - 1) - 2^(iV-l)) 



where the sum is over all permutations ir and ir' which differ by a single transposition involving 
7r(0). Combining this with Lemma 1 we obtain 



I6h* 



1 



^2 



Ah* 



which is what we wanted to prové. 



□ 



What we have just shown is that no continuous time algorithm of the form (|14| can find the 
minimum of -£/p )7r with a constant success probability for even a fraction eiV! of all permutations 
7r if T is o{\fN). A typical permutation ir yields an Hp n with no structure relevant to any fixed 
and the algorithm can not find the ground state of iíp j7r efficiently. 

To illustrate the nature of this failure for the quantum adiabatic algorithm for a typical per- 
mutation, consider again the decoupled n bit problem with h(z) given by and Hb given by 
(fT2l . The lowest curve in FIG. 2 shows the ground state energy divided by n as a function of t. 
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FIG. 2: The scaled ground state energy E/n for a quantum adiabatic algorithm Hamiltonian of a decou- 
pled problcm. The lowest curve corresponds to the original decoupled problcm. The upper "triangular" 
curves correspond to singlo instances of the rt-bit decoupled problem, where the problcm Hamiltonian was 
scrambled. 



(Since the system is decoupled this is actually the ground state energy of a single qubit.) We then 
consider the n bit scrambled problem for different vàlues of n. At each n we pick a single random 
permutation ir of 0, . . . , (2 n — 1) and apply it to obtain a cost function h{iï~ l {z)) while keeping 
Hb fixed. The ground state energy divided by n is now plotted for n = 9, 12, 15 and 18. From 
these scrambled problems it is clear that if we let n get large the typical curves will approach a 
triangle with a discontinuous first derivative at t = T/2. For large n, the ground state changes 
dramatically as t passes through T/2. In order to keep the quantum system in the ground state 
we need to go very slowly near t = T/2 and this results in a long required run time. 



IV. CONCLUSIONS 

In this paper we have two main results about the performance of the quantum adiabatic al- 
gorithm when used to find the minimum of a classical cost function h(z) with z = 0, . . . , N — 1. 
Theorem 1 says that for any cost function h(z), if the beginning Hamiltonian is a one dimensional 
projector onto the uniform superposition of all the \z) basis states, the algorithm will not find the 
minimum of h if T is less then of order y/N. This is true regardless of how simple it is to classically 
find the minimum of h{z). 

In Theorem 2 we start with any beginning Hamiltonian and classical cost function h. Replacing 
h(z) by a scrambled version, i.e. fcH( z ) = h(n(z)) with 7r a permutation of to N — 1, will make it 
impossible for the algorithm to find the minimum of hfà in time less than order y/N for a typical 
permutation n. For example suppose we have a cost function h(z) and have chosen Hb so that the 
quantum algorithm finds the minimum in time of order log N. Still scrambling the cost function 
results in algorithmic failure. 

These results do not imply anything about the more interesting case where Hb and Hp are 
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structured, i.e., sums of terms each operating only on several qubits. 
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APPENDIX A: TRANSITIONS IN A TWO LEVEL SYSTEM 

Let us consider a two level system with Hamiltonian 

H( S ) = E ( s )\M s ))(Ms)\ + EiOOl&OOX&OOI, 

which varies smoothly with s = t/T. Here \(j>o{s)) and \(pi(s)} are orthonormal for all s. The 
Schròdinger equation reads 



i—W)=TH(s)W). 



The two energy levels in the system are separated by a gap 

g (s)=E 1 (s)-E (s), 

which we assume is always larger than 0. Let us introduce 9 (with the dimension of energy) as 



9( S ) 



g(s') ds', 



and let 



(Al) 



We pick the phases of \<pi{s)) and \</>o(s)) such that (4>i(s)\£\(f)i{s)) = (</>o(s)|^|^o(s)) = 0. 
Plugging (Al) into the Schròdinger equation gives 



dcp 
ds 
dci 

~ds~ 



ds 
d 

ds" 



0f 



or equivalently, 



where 



d 

de 



dço _ iT8 f * 

de ~ Cie 1 ' 

de ° h 



dH 
~dJ 



dH 
ds~ 
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Now let 9(1) = 8. We have c\(0) = and we want the transition amplitude at s = 1 which is 



coé Te fd9 



-cof 



1 

T 



JT9 

ic fe 



+ 



o 

iTd 



ÍT 



ÍT9 /' dc , , df\ 
6 'd^ / + CO d^ ld0 



are bounded so we 



Now |co| < 1 and |ci| < 1. As long as the gap does not vanish and 

have that |ci(#)| = O )• The probability of transition to the excited state for a two-level system 
with a nonzero gap is thus 



\ Cl (e)\=o(T- 
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